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Abstract 



The electromagnetic form factors of hadrons at large momentum transfer have been the subject of intense 
theoretical and experimental scrutiny over the past two decades, yet there is still not a universally-accepted 
framework for their description. This review is a synopsis of their current status at large momentum transfer. 
The basic theoretical approaches to form factors at large momentum transfer are developed, emphasizing 
the valence quark and Feynman (soft) pictures. The discussion includes the relation of these descriptions 
to the parton model, as well as the roles of factorization, evolution, Sudakov resummation and QCD sum 
rules. This is followed by a discussion of the experimental status of pion and nucleon elastic form factors and 
resonance production amplitudes in the light of recent data, highlighting the successes and shortcomings of 
various theoretical proposals. 
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3 Experimental Status of Hadronic Form Factors 



1 Introduction 

1.1 Hadronic Form Factors 

Exclusive electromagnetic form factors are a source of information about the internal structure of hadrons. The 
coupling of an elementary particle to the photon is determined by only a few dimcnsionlcss parameters, for 
example its total charge and magnetic moment. For a composite particle, however, these constant coefficients 
are replaced by momentum-dependent functions, the form factors, which reflect the distribution of charge 
and current, and hence the internal structure of the particle. A familiar analysis in nonrelativistic quantum 
mechanics relates the electromagnetic form factor directly to the Fourier transform of the the charge density. 
Rclativistic behavior also depends very much on the nature of the hadronic state. 

High momentum transfer suggests high resolution, so hard elastic scattering is a natural way to study 
the detailed internal structure of hadrons. Experiments in elastic electron-proton scattering showed long ago 
the famous dipole behavior of the nuclcon electromagnetic form factors in terms of momentum transfer Q, 
1/(1 + Q 2 /M 2 ) 2 , with M ~ 0.71 GeV 2 §. 

Since then, many subsequent experiments studied this and related reactions. Their influence on our un- 
derstanding of the strong interactions themselves, however, has been somewhat overshadowed by that of the 
high energy inclusive reactions. The discovery of approximate scaling in deeply inelastic scattering, and its 
explanation in terms of the parton model, opened a more direct and efficient avenue to study the quarks them- 
selves, since inclusive rates decay much more slowly with momentum transfer. Nevertheless, form factors at 
large momentum transfer remain an important window to quark binding in hadrons. 

In this review, we will concentrate on electromagnetic form factors and resonance production amplitudes, at 
large momentum transfer, in the light of perturbative quantum chromodynamics (QCD). QCD itself has enjoyed 
so many successes, and explains so many and varied experimental results, that it is universally recognized 
as "the" theory of the strong interactions. Yet, the single most basic fact of the theory, the binding and 
confinement of the elementary degrees of freedom, the quarks and gluons, into hadrons, is still not described 
in detail. Because of the property of asymptotic freedom at short distances, perturbative methods must be 
relevant in some degree to elastic scattering at large momentum transfer. Because the binding of hadrons is a 
long-distance effect, nonperturbative effects must play a crucial role as well. The description of electromagnetic 
form factors requires the consistent analysis of both length scales in a single process. This, and the light that 
will be shed on hadronic structure by a truly successful treatment of this problem makes the study of form 
factors attractive. We note that electromagnetic form factors are part of the large class of exclusive hadronic 
amplitudes, which also describe, for example, both proton-proton elastic scattering and the exclusive decays of 
heavy mesons. Although many of the methods developed below have wide applications in this larger class, we 
decided to restrict our discussion to form factors, in the hope of improving its focus. 
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In the remainder of this section, we discuss what we can learn from reasoning based on the parton model. 
Here, and in most of the following, we assume very high momentum transfer, so high that parton masses may 
for the most part be neglected. We use parton model insights to identify quark counting rules, and as an 
inspiration for factorization of long- and short-distance effects in exclusive processes in terms of wave functions. 
In this section, we give primarily intuitive arguments, and concentrate for simplicity on the pion. In Sec. ^, 
we discuss some of the central results of the QCD treatment of form factors, including the evolution of wave 
functions, the behavior of the asymptotic pion form factor, and QCD sum rules for moments of wave functions. 
These topics are somewhat more mathematical, but we have attempted to motivate technical arguments with 
physical intuition. We close this section with a brief summary of results relevant to baryons, especially helicity 
conservation, the derivation of form factors directly from QCD sum rules, and a few phenomenological models 
for moderate-Q 2 behavior. We go on to review the central experimental results for pion, nucleon and resonance 
production form factors, to assess the successes and failures of QCD treatments of elastic scattering, and to 
explain the controversies that have enlivened this active field of inquiry. To anticipate, we will see that the 
current state of the data is not adequate to resolve the primary theoretical controversies. 

1.2 Partons and Factorization 

1.2.1 Partons. The perturbative treatment of hard exclusive processes assumes a partonic description of the 
participating hadrons. The general discussion is closely related to the parton model of inclusive processes 
such as deep-inelastic scattering. The celebrated premise of the parton model, justified and systematically 
extended in QCD, is that inclusive processes are determined by the distributions fi/h(x), which are the proba- 
bilities for pointlike, constituent partons i to carry fraction x of the momentum of hadron h, summed over all 
other partonic degrees of freedom. An exclusive form factor, on the other hand, reflects the coherent scattering 
of a hadron by an electroweak current. Even at large momentum transfer, it may depend on states of definite 
partonic content. In fact, at high enough energies, exclusive amplitudes are dominated by hadronic states with 
"valence" quark content, qq for mesons and qqq for baryons. This is despite the fact that, in its own rest frame, 
each hadron is a complicated, ever-shifting superposition of partonic states. Let us discuss first how such a 
partonic picture of hard exclusive scattering emerges. 

At sufficiently high momentum transfers in either hadron-lepton or hadron-hadron scattering, the relative 
velocities of all participating particles are nearly lightlikc. Under this condition, the quantum processes that 
bind the constituents of a hadron are highly time-dilated in the rest frames of the remaining particles, both 
incoming and outgoing. Correspondingly, time dilation lengthens the lifetime of these states, and "freezes" the 
partonic content of this hadron as "seen" by the other particles. Also, as relative velocities approach the speed 
of light, the time during which the hadrons remain in contact, and during which momentum can be transferred, 
decreases. In fact, we can always find a frame in which any pair of particles are in contact for a time that 



between long-distance, hadronic binding and short-distance momentum transfer. This incoherence between soft 
and hard physics implies that we may consider each hadron to consist of a definite partonic state during the 
entire collision process. This picture is illustrated for electron-pion scattering in Fig. |l|a, in which long-time 
dynamics, described by a distribution of valence quarks 4>i n , produces a "valence" quark-antiquark state. The 
distribution <f> is often referred to as a "wave function" . The partons of this state in turn exchange momentum 
with an electron in a short-distance process T . At a later time, they reform a pion, through wave function (j) ou t- 




Under these conditions, we expect a lack of quantum interference 



0; 




(a) 



(b) 



Figure 1: Electron-pion elastic scattering, a) Valence PQCD picture, b) Feynman mechanism. 



1.2.2 Factorization. We summarize the above considerations for an arbitrary exclusive amplitude M by a 
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schematic expression in which short-distance momentum transfer is factorized from the long-distance hadronic 
binding, 

M(pi-pj) =Y[(f) outtj (nj) ®T(nj,m) ®\\(j)in,i{ n i) • (!) 

3 i 

Here, the labels i and j refer to hadrons in the incoming and outgoing states, respectively. (j>(n) is the wave 
function that describes the amplitude for a pion to be found in partonic state n, and T(rij, n,) is a perturbative 
function that describes the hard scattering between the partons (and leptons). The symbol ® indicates a 
convolution, that is, a sum or integral over the parton degrees of freedom that correspond to states rii and rij . 

A factorized expression like Eq. (|l|) has two fundamental properties. First, the nonperturbative wave func- 
tions are universal within a class of exclusive amplitudes. This connects otherwise disparate processes, such as 
the pion electromagnetic form factor and pion-pion elastic scattering || . Second, the factorization of long- from 
short-distance dynamics implies consistency conditions that enable us to compute the amplitude's dependence 
on the momentum transfer. These are usually referred to as "evolution" equations, examples of which we shall 
discuss below. The details of the convolution ®, and the derivations of evolution equations, depend on the 
process in question, but one example will suffice to motivate Eq. (Q) and to illustrate the range of possibilities, 
the electromagnetic form factor of the charged pion. We will review the classic perturbative QCD analysis of 
this form factor j|, ||, ^| , and also introduce a treatment of its "Sudakov" effects jj], |s| , whose importance will 
become clear below. 

1.2.3 Valence PQCD and the Feynman Mechanism. The convolutions ® in Eq. ([!]) in principle include 
sums over states with arbitrary numbers of partons. As indicated above, however, at very large momentum 
transfer, the valence state, with the fewest partons, dominates. We shall refer below to its contribution as 
"valence perturbative QCD" (valence PQCD). This is a somewhat unconventional usage; indeed, what we call 
valence PQCD is more commonly referred to simply as "PQCD". But this approximation does not exhaust the 
use of perturbative methods in form factors at large momentum transfer, and to call it simply "PQCD" is a 
little misleading. 

There are, of course, many contributions from states with more than the valence partons. For the most part, 
they are expected to decay rapidly with increasing momentum transfer, relative to the valence states. There is 
an exception, however, corresponding to states in which one parton carries nearly all of the hadron's momentum, 
while all other partons are soft. It is plausible that such a state could contribute to elastic scattering, because 
all of its partons except for one have long wavelengths. They may then overlap strongly with wave functions 
moving in any direction. When the single, hard parton scatters elastically, the soft partons from an incoming 
hadron may combine with the outgoing hard parton to form an outgoing hadron. This is illustrated for the pion 
electromagnetic form factor in Fig. |l|b. It is known as the "soft" or "Feynman" mechanism for elastic scattering. 
Nevertheless, the Feynman mechanism contains a hard scattering, which may, in principle, be factored from the 
interactions of soft partons, and treated with the methods of PQCD. For instance, in || [To) it was analyzed for 
pion and nuclcon form factors. This PQCD investigation, unfortunately, has not yet been developed extensively 
in the literature, and although it seems clear that the soft mechanism does not contribute at asymptotically 
high momentum transfer, at what scale it becomes negligible is not well understood. We shall come back to 
the role of the soft mechanism often below, however, because its contribution may be studied directly in the 
valence state, using Sudakov resummation, and in QCD sum rules, and, indirectly, in models of nonperturbative 
hadronic structure. 

1.2.4 The Pion Form Factor and Quark Counting. The electromagnetic form factor of a pion is 
specified by 

( P 2+Pi)^F„(Q 2 ) = (n(p 2 )\ M0)\*M . (2) 

where = X)/ e /9/7^9/ ^ s the electromagnetic current, expressed in terms of quark fields g/ of flavor / and 
electromagnetic charges e/. We neglect particle masses, and examine this process in a "brick- wall" frame, 
in which pi is in the plus 3 direction, and recoils as pi in the minus 3 direction under the influence of the 
electromagnetic current J. Such a momentum configuration is most naturally described in terms of light-cone 
variables, which for any vector are w ± = 2~ 1 / 2 (i>° ± v 3 ). In these terms we have 

p+ = Q/y/2, p{ = 0, p 2 = Q/y/2, p+ = . (3) 
The overall momentum transfer is (p 2 — Pi) 2 — ~2p^p 2 = —Q 2 . 
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The valence PQCD portrait of this process is shown in Fig. ||. Fig. ||a represents the pion in its valence 
state, consisting of a quark and an antiquark. The variable x denotes the fraction of the pion's momentum 
carried by the quark and 1 — x by the antiquark. In the chosen frame, we expect the pion to be Lorentz 
contracted in the direction of motion, as shown, so that the pair is localized in this direction. On the other 
hand, we expect the partons in any virtual state to be more-or-less randomly distributed in the transverse 
extent of the pion's wave function, since the boost from the rest frame to the frame under consideration leaves 
transverse positions unchanged. This will have important consequences below. Similarly, the off-shellness, and 
the transverse momenta of the pair in Fig. ^| are boost-invariant, and we take these quantities to be fixed, and 
negligible compared to both xp% and (1 — x)p%. Correspondingly, the transverse components of their velocities 
vanish as Q — > oo, and we neglect them as well. It is necessary that 1 > x > 0, so that both partons travel in 
the same direction as the hadron that they represent. Fig. ^a also shows an incoming, off-shell photon, carrying 
momentum q. 



1/QJ 

1 



A 



xp , 



(1-x)p 



Jt(p,) Jt(p 2 ) 



V 



A 




(a) 



(b) 



Figure 2: Pion electromagnetic form factor in valence PQCD. 



Fig. |2jb shows the state of the system after the action of the current that absorbs the photon, in which the 
pair moves in the opposite direction. Eventually, the pair will fill out the full spatial extent of the pion, which 
is again Lorentz contracted. To form the pion, however, their momenta must be parallel, and each must carry 
a positive fraction of P2, as shown. 

An alternative picture relies on the "infinite momentum frame" (IMF), in which all participating particles 
move in the same direction, with energies Ei 3> Q. In this frame, all momentum transfers are transverse. 
Its main attraction lies in the conjecture that quantization formulated in an IMF simplifies the treatment of 
confinement in QCD [pd| . 

In the process depicted in Fig. ^, the quark undergoes a momentum transfer xyQ 2 , and the antiquark 
(1 — x)(l ~ y)Q 2 , with x (y) the fractional momentum of the quark in the incoming (outgoing) pion. This must 
take place during the time that the wave functions of the incoming and outgoing pions overlap, that is, on a 
time scale that vanishes as l/Q. The uncertainty principle requires that both members of the pair must be 
localized within l/Q of each other and of the action of the current, as indicated in Fig. ||a. This restriction 
shows, first of all, that not all details of the valence state wave function are relevant to exclusive scattering. We 
do not need the full two-particle state; we only need the probability for the members of the pair to be within a 
transverse distance of l/Q of each other. We shall assume that this probability is simply a function of x times 
the geometrical factor l/Q 2 . This "scaling" of the wave function in x is not exact in QCD; we will compute 



corrections to it when we discuss evolution below in Sec. 2.2 



Along with our assumption of incoherence, scaling enables us to estimate the Q-dependence of the form 
factor. For, if long- and short-distance processes are incoherent, the cross section for elastic scattering of a 
pion is essentially the product of the cross section for the elastic scattering of a point-like scalar particle, times 
the probability for internal processes to produce a virtual state in which both partons in the valence state are 
within l/Q of each other in transverse distance. Thus, we have 



Cel, point X F 2 (Q 2 ) ~ (Tel, poin t X (l/Q 2 ) 2 



so that 



FAQ) ~ (i/Q 2 



(4) 



(5) 



Results of this sort, based on incoherence, scaling and geometrical estimates, are known as "quark counting" 
[ |l2] , fl3|| . Quark counting rules give for an arbitrary exclusive process involving nh hadrons, 



a{Q< 



)had — Cr(Q 2 ) po i n t 



(Q 2 



7, 



(G) 
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where n q is the total number of quarks and antiquarks taking part in the process, and / depends on dimcnsionless 
variables. 

iFrom Eq. (§), we see that interactions involving more than the minimum number of partons - say, a gluon 
in addition to the pair - are suppressed by a power of Q, because as Q grows, the likelyhood of finding more 
than the minimum number of particles within 1/Q of each other falls as Q~ 2 for each additional particle. 

We note, however, that in the limits x,y — > 0, 1, our process describes the elastic scattering of an on-shell 
quark (antiquark) with nearly all of the pion's momentum. The remaining, soft antiquark (quark) has long 
wavelength, which overlaps with both the incoming and outgoing wave functions. This is the intersection of 
valence PQCD with the Feynman mechanism. 

In the next section, we shall turn to the field-theoretic treatment of the pion's form factor, and shall see 
how these features of the parton model are realized within QCD. 



2 Form Factors in QCD 

2.1 The Factorized Pion Form Factor 

2.1.1 Convolution in Fractional Momenta. We are now ready to turn to the pion form factor in valence 
PQCD j|, ||, |(| . The parton model discussion of the previous section suggests that the pion form factor can be 
written, following Eq. (|l|), as a sum over wave functions involving only quark momentum fraction. We denote 
these as x and y for the incoming and outgoing pions, respectively. We then have the following representation 
for the form factor, 



Fir(Q 2 ) = / dx dy ^(y, /i 2 ) T(y,a;,Q 2 ,^ 2 ) (j} n (x,fi 2 ) . 
Jo 



(7) 



Here, ^(x, [i 2 ) is the valence-state wave function describing a quark with fraction x of the pion's momentum. 
T(x,y,Q 2 , /j, 2 ) describes the hard scattering of partons. It is a perturbative expansion in the strong coupling 
at scale /x 2 (a s (/i 2 )), and is free of infrared divergences order- by-order in perturbation theory. At lowest 
order, it is given by the diagrams shown in Fig. |3|a. Since the incoming and outgoing pairs are each at a 
tiny transverse separation, orbital angular momenta are negligible, and partonic helicities must sum to zero. 
The pairs of incoming and outgoing external lines in the diagrams are thus projected onto Dirac matrices that 
represent these helicity-zero pairs. [|| At the same time, because masses are neglected, helicities are conserved 
in perturbation theory, and hence in the hard scattering, to all orders. This has important consequences for 
hadrons with spin. An exercise in dimensional counting shows that T has dimensions (mass) -2 , and hence 
scales as 1/Q 2 . Its explicit lowest-order form is 



T H = l6-KC F a s (fi 2 ) 
where Cf 



1 



1 



1 



3xyQ 2 3 (1 - x)(l - y)Q 2 
(N 2 - l)/2N = 4/3 in QCD. 



(8) 
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(a) 



(b) 



Figure 3: Graphical contributions to T: (a) Lowest order; (b) Typical one-loop correction. 

In perturbation theory, it is possible to show that the valence PQCD result Eq. (0) is a theorem, which 
describes the behavior of F^iQ 2 ) at large Q 2 to all orders in a s . Corrections are suppressed by powers of Q 2 , 
including those due to the Feynman (soft) mechanism described above JE], |(], 0. Helicity conservation in T is 
also valid up to similar corrections. Space does not allow a discussion here of technical aspects of this proof or 
of the calculation of T beyond lowest order. The basic technique is already illustrated by a typical one-loop 
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correction, Fig. |b. (For the pion form factor, the full one-loop calculation has been performed explicitly 
[ |l4"[ We keep only those contributions to T from Fig. ||b where all quark and gluon lines are off-shell by 
at least the renormalization scale, n 2 . T then depends on only two momentum scales, Q and /i. Alternately, 
we can think of /x 2 as the minimum transverse momentum carried by lines in T. Suppose, now, that we choose 
fi = Q. By the uncertainty principle, this corresponds to including in T only lines that are within 1/Q of each 
other in transverse distance, as we anticipated in our discussion of Fig. || above. By choosing Q — /x we get the 
extra benefit of expanding T in terms of the small parameter a s (Q 2 ). Thus, /x = Q will be our default choice 
of scale below, although other choices may sometimes offer special advantages. 

2.1.2 Transverse Degrees of Freedom in F„. Our next exercise in factorization is to return to Eq. ([?]), 
taking into account transverse degrees of freedom. Remaining in the valence picture, we recall that the pair in 
the incoming and outgoing pions are not literally at a point, but are separated by transverse vectors bi when 
they undergo the hard scattering, where i = 1 (2) for the incoming (outgoing) pion. Again, Q 2 = (p2 — Pi) 2 - 
The wave functions in Eq. (Q), which we now denote V, are characterized by both fractional momenta and 
transverse separation, and the form factor is reexpressed as a convolution in both Q , 

F^(Q 2 ) = J dxdy J V(y,b 2 ,P2,fJ>) 

x T(y,x,pi,b,n)V{x,bi,Pi,IJ) > (9) 

where T is a new hard-scattering function. Beause we integrate over the variables bi conjugate to transverse 
momenta, fj, does not play the role of a "transverse momentum cutoff", as in Eq. (Q), but is simply the renor- 
malization scale. On the other hand, the wave functions V depend upon the momenta pi and they, along with 
T, are not individually Lorentz invariant. The requirement of Lorentz invariance in the complete amplitude F v 
will lead to evolution equations below. We emphasize that, summed to all orders, Eq. (9) is equivalent to Eq. 
(f7|) at leading power in Q 2 . Depending on the details of V, however, it differs from Eq. (|7) in nonleading powers 
of Q 2 in general. 

We next explore the relation between the two factorization procedures a little further. Intuitively, we expect 
that the wave function V near 6 = 0, that is at small separation for the pair, is related to the distribution 
amplitude (f>^. Specifically, it is not difficult to show that |Q 

V(x,b=l/(j,,pi,n) ~ 4>(x,fi 2 ), (10) 

up to corrections that are suppressed by the strong coupling evaluated at the factorization scale fi. The Lorentz 
noninvariance of the V disappears in this limit. 

Let us now compare Eq. (||) to the classic expression, Eq. If Q 2 is large enough, we expect, according 
to our discussion above, that T in Eq. @ is concentrated near b ~ 1/Q, so that by Eq. ([n]) V may be replaced 
by (/>„. In this limit, the two expressions are equivalent. A closer look at T in Eq. (0), however, shows that it 
actually corresponds to a localization in transverse space only at the scale (xyQ 2 ) -1 . When x or y vanishes, 
the hard scatteri ng " spreads out" in transverse space, and violates the original assumptions of the partonic 
discussion of Sec. |l.2| .l above, and the reaction is defined by the Feynman mechanism. Also note that if b is 
large, the neglect of orbital contributions to hclicity is no longer justified, even if hclicity is conserved in the hard 
scattering [[l6| [l?], |l^|. The contribution of the "end-point regions" x, y — > 0, 1 (the equivalent of the Feynman 
mechanism for qq states) depends on the details of the wave functions 0, but it poses a problem, unless Q is 
very large (II], |o| We shall see shortly that the use of the modified factorization in Eq. (||) serves to stabilize 
the valence PQCD picture of scattering at somewhat lower Q 2 than in Eq. ^ ||. To see how this comes about, 
we turn now to a discussion of evolution, as derived from the factorization formulas Eq. (Q) and Eq. (||). 

2.2 Evolution and Asymptotic Behavior 

Eqs. (Q) and (^) for the elastic form factor are both convolutions of functions that depend upon arbitrary 
choices: the renormalization scale /x in the former case, and the Lorentz frame in the latter. In fact, a great deal 
can be learned from these parameters, through their role in the factorization formulas. Among other things, it 
will allow us, in the following subsection, to give an explicit expression for the asymptotic behavior of <f> n and 
the form factor at high momentum transfer. 

2.2.1 Evolution. Consider Eq. (0) for F^{Q 2 ). The physical form factor, of course, cannot depend upon 

fx: 

^F 7r {Q 2 )=i). (11) 
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Equivalcntly, in terms of the hard-scattering and wave functions, 



= / dxdy 



d T4>„[x) + (p- K {y)—(p, r {x) 



<PAy) T - 



dfi 



(12) 



This expression may be treated by separation-of- variable techniques. d(f) n (y, fx )/dfX, for instance, may depend 
upon the variables y and p, 2 only, the latter only through a s (p 2 ) (since there are no other dimensionless variables 
available.) In fact, its derivative with respect to p 2 must be perturbatively calculable, because changes in p 



shift contributions from lines that are off-shell by order p 2 between <\>^ and T. (See Sec. 2.1.1).) The most 
general form that satisfies these requirements is itself a convolution Ifj) : 



dp 



The kernel V is a distribution, rather than a simple function of y and z, but its integral with any smooth function 
is finite. Given the convolution form Eq. (]?]) for the form factor, the evolution equation (|l^) holds to all orders 
in a s (p 2 ). Its explicit one-loop form is simply the coefficient of InQ 2 in the sum of one-loop corrections to the 
hard scattering, such as Fig. The kernel V is known up to two loops (2l). We shall not exhibit its explicit 
form, but only note that, with the one-loop V, Eq. (O) may be solved explicitly. The most general solution is 
an expansion in Gegenbauer polynomials C n 3 / 2 || J22| . 

<Mx, M 2 ) =x(l-x) y £ a n Cl' 2 {2x l)(ln , (14) 

n>0 

with /?2 = (33 — 2n/)/12 the one-loop coefficient of the QCD beta function, the 7„ known anomalous dimensions 
and the a n arbitrary coefficients. 

Space allows us to make only a few observations on this fascinating result: (i) The a n are linear combinations 



of matrix elements, identified in Sec. 2J3 below; (ii) 70 = 0. This is because the n — wave function, 4>o(z) = 
aoz(l — z) gives zero when integrated with the one-loop kernel in Eq. (|l3|). We shall refer to this "asymptotic" 
form of the pion wave function many times below; (iii) for n > 0, all j n > 0, which implies that as p 2 — > 00, all 
^-dependence in (|l4|) that is not in the form of the asymptotic wave function decays, albeit only logarithmically. 

2.2.2 Sudukov Resummation. Turning now to factorization in transverse space, we see that the fac- 
torization Eq. (^) suggests another evolution equation, this time in the momentum scale Q, which enters the 
wave functions through (non- invariant dependence on) the momentum vectors pi. This equation will allow us to 
resum perturbative logarithms of the form ln(Q6), with b the distance between the hard scatterings, an example 
of "Sudakov resummation". The derivation of this equation is given in a related context in H. Here, we shall 
content ourselves with a physical explanation and the basic results. In brief, the effect of the resummation will 
be to suppress the nonperturbative contribution to F n , and thus to extend valence PQCD to lower Q 2 . 

For large momentum transfer, the dynamics of elastic scattering strongly disfavors configurations in which b is 
large. The physical reason for this result is that an isolated accelerated charge must radiate, by correspondence to 
classical gauge theory. As b grows, the two charges associated with quark and antiquark become more isolated, 
and have correspondingly more tendency to radiate gluons. In elastic scattering, however, such radiation is 
forbidden by definition. Perturbatively, this manifests itself in the presence of double- logarithmic ("Sudakov") 
corrections of the form a™ (/z) ln 2n (bQ). We therefore expect that the double logarithms at large b will suppress 
configurations for which the charges are separated far enough to couple strongly to radiated gluons. Because 
the effect is essentially classical, it is necessary to sum to all orders (take the limit of large quantum numbers) 
to make this suppression manifest. 

In this case, an evolution equation is derived from the independence of expressions like Eq. (|^) of the choice 
of inertial frame. An infinitesimal Lorentz transformation changes the arguments of the P's and of T, but 
otherwise leaves the amplitude invariant. A full derivation (sec Q|; the reasoning there is an application of a 
method first developed in Ref. |23j]) requires more analysis of b and Q dependence in T than we have room for 



here. The result is the following evolution equation, which takes the place of Eq. (13). Taking p + = Q in the 
center-of-mass frame, we have, 

Q-^P{x 7 b,p, n) = \K(bp) + G(x, Q//x)] P(x, b,p, p) , (15) 
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in which the functions K and G may be computed in perturbation theory. K depends only on the "infrared" 
variable b, and G on the "ultraviolet" variable Q. 

The details of the solution to this equation is straightforward, and may be found in H. The result is striking: 



, b;p,fi) = e Si*AQ,rt ^( s> 1/6 3) + ( a 2 s (l/b)) ) , (16) 

where <p n is the usual light-cone wave function for the pion, now evaluated at fi — 1/b. The Sudakov exponent 
S strongly suppresses the wave function at large b, through the summation of double logarithms of bQ per loop, 

f xQ dn' a s (fx') fxQ\ 

S = Cf — ; In + x «-» 1 — x + . . . , (17) 

Ji/b M V M / 

where we have suppressed terms with fewer logarithms per loop. Note in particular that within the integral 
the perturbative coupling runs with the variable //, so that the Sudakov exponent S diverges at b = 1/Aqcd- 
When Q ^> Aqcd, the exponent is large, and the suppression great, whenever b 3> l/Q, even for b <C 1/Aqcd- 
The quark- antiquark state with opposite helicties again dominates in this limit. The suppression of large-6 con- 
figurations has many applications to hadron-hadron reactions, and helps justify the concept of "transparency" 
in hadron-nucleus scattering |24| ]. 

2.2.3 The Asymptotic Form Factor. We are now ready to discuss one of the central results of the 
perturbative treatment, the asymptotic behavior of the pion electromagnetic form fa ctor . We begin by recalling 
that the natural choice of scale in the factorized expression Eq. (Q) is Q = /x (See Sec. 2T.1). For Q large enough, 



then, the wave function will be dominated by the ciq term in its expansion ( |l4|) , and T will be well-approximated 
by its lowest-order contribution, Eq. The x and y integrals in (Q) are then simple, and the only remaining 
uncertainty is in a factor of a 2 . 

To fix do, we observe that the decay of the charged pion though the weak interactions may be treated by 
the same method of factorizing hard and soft degrees of freedom. In this case, the hard interaction is at a scale 
of the order of the W-mass, and the wave function of Eq. is again dominated entirely by its qq coefficient. 
Then, defining the pion decay constant by (with f„ ~ 93 MeV) 

<0|d(0) 7 M (l - 7fe)tt(0)|7r(p)) = -vVA , (18) 

we may identify ao = VSf^ in Eq. (pi). (More generally, we have 3f w /y/N with N the number of colors). This 
result, along with properties of the anomalous dimensions ("y„ > for n > 0), allows us to identify the large /i 2 
(or Q 2 ) behavior of the pion's quark wave function: 

<Mx,/i 2 ) -» y/3f«x(l-x). (19) 

This is generally referred to as the asymptotic wave function of the pion. We emphasize that it is model- 
independent. 

Substituting Eq. ( fl9|) into Eq. (0), and using the lowest order hard-scattering function T (Eq. @) with 
/i = Q, we find an elegant expression for the pion form factor at high energy, which is valid up to corrections in 
a.(Q)~l/ln(Q) |,|, 

12/2ttCWQ 2 ) 
Q 2 

2.2.4 Sudakov Resummation for F n . With an eye to contributions for which the pair is widely separated, 
we may also use the Sudakov-resummed transverse wave function ([ll]) in Eq. (0) , to get 

FAQ 2 ) - f o dxdy J -0^ MV, V& 2 ) e-^ 6 >^) 

x T(y, x, Q, b, ^) e -s(*,*,Q,i>) ^( Xj i/ftS) , (21) 

where we have simplified to a single transverse separation ||, ^5|. The Sudakov exponential in this expressions 
factor suppresses contributions from b 1/Q. The natural scale of the coupling in T is fi ~ 1/6, even in the 
end-point region. Perturbation theory thus remains self-consistent, by the dynamical suppression of the overlap 
region of valence PQCD and the soft mechanism. For moderate Q 2 , however, F n still receives substantial 
contributions from relatively large b < 1/Aqcd- In this region, ( pl| ) should be thought of as a valence PQCD 
model for F„ ■ Form factors computed according to each of these procedures will be confronted with the data in 
Sec. H below. 



F*m= '* X?'™ ' ■ (20) 
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2.3 Wave Functions and Nonperturbative Analysis 

In the factorized picture of elastic scattering we treat hadrons as superpositions of states, each with definite 
numbers and positions (or momenta) of partons. Also, as we have seen, it is the states with the fewest partons 
- the "valence states" that dominate exclusive processes at sufficiently high Q 2 . Relativistic valence wave 
functions for valence states may be identified with matrix elements that connect single-particle states of defi- 
nite hadron momentum \h(p)}, with the hadronic vacuum |0) by the action of fields that absorb the relevant 
valence quanta. The analysis of these matrix elements can lead to valuable nonperturbative information, which 
supplements the purely perturbative results outlined above. 

2.3.1 Matrix Elements. The light-cone wave function in position space for the valence state of a ir + may 
be defined in terms of the matrix element of an up quark field u with a conjugate down quark field, 

*(* -p,z 2 ) = (0\d{0) 7 + 7 5 w(z)|^+(p)) , (22) 

where p is taken in the plus direction. In the following, we shall generally neglect hadronic masses. We recall 
that 7 ± = l/y/2 (7 ± 7 3 ). So that ^(z ■ p,z 2 ) may have a natural interpretation in terms of independent 
measurements of the up and antidown quark fields, we choose the separation between the two fields to be 
spacelike, z 2 < 0. The Dirac structure 7 + 75 projects out precisely the zero-helicity combinations of the quark 
and antiquark fields. 

As defined, the wave function \P of Eq. (|2^) is gauge-dependent. A common choice of gauge for the gluon 
field A is A + = for a pion moving in the plus direction. Alternately, we may connect the fields d(0) and u(z) 

by a path-ordered exponential in the direction Pexp J dtz^ A^(zt) , with A^ expressed as a matrix in the 

quark representation. 

The momenta of the partons in the valence state may be fixed by taking Fourier transforms. For (j) v (x), we 
fix the fractional momentum of the quark to be xp in the pion's direction of motion, and integrate freely over all 
of its other components (and hence those of the antiquark) by setting z + = zt = 0, and taking the transform 
of ^> with respect to z~ , 

/OO 7 — 
^_e fa ^ + *(z-p+,0) A+=0 . (23) 

Here 4>n( x ) depends explicitly on the renormalization scale /1 2 , because the limit z + , zt — > 0, which takes z^ to 
the light cone, z 2 = 0, is singular. Defined in this fashion, ^(x) is referred to as a light cone wave function. 

Readers familiar with the QCD analysis of deeply inelastic scattering (DIS) will recognize a similarity 
between the valence quark wave function given by Eqs. ( p2| ) and ( |23] ) and the inclusive parton distribution 
density in a hadron. Note, however, that while a parton distribution in DIS is a probability, (j)^ is an amplitude. 
Thus, although we know the behavior of our light-cone wave functions at very large [i, they might evolve slowly 
to this form, and we would like further information on their properties for intermediate values of /i. A direct 
approach is to compute the relevant matrix elements using the methods of lattice QCD. Moments of proton 
wave functions have been computed in this fashion (2^, ^7|, and more work may be antitipated in the future. 
Direct, nonperturbative information on the wave functions may also be found using instanton models of the 
QCD vacuum ||. 

The traditional approach to derive extra, nonperturbative knowledge on wave functions has been the use of 
QCD sum rules to determine their moments of with respect to x. We shall discuss light-cone wave functions 
only, but we note that sum rules have recently been applied to wave functions with transverse degrees of freedom 
© 0. 

2.3.2 Sum Rules for Wave Functions. QCD sum rules |29] have many applications, whenever a non- 
perturbative quantity can be related by analyticity to the integral of a Green function (vacuum expectation 
value of a time-ordered product of local fields) over a range of highly virtual momenta. When this is the case, 
perturbation theory, supplemented by the operator product expansion (OPE), may be used to calculate the 
integral of the Green function, from which the value of the matrix element may then be inferred. 

In the following, we show how QCD sum rules may be used to obtain the moments of wave functions, 
parameterized in terms of experimentally fitted gluon and quark vacuum condensates [ ^9| . Using this technique, 
Ref. I32J obtained the following simple wave function, 

0cz (x, Mo) = 5V3/ T a;(l - x){\ - 2x) 2 , (24) 

with fio ~ 0.5 GeV. This result became a common test case for many subsequent authors. This wave function 
is plotted, along with the asymptotic wave function, in Fig. W. Compared to the asymptotic expression of Eq. 
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|19|), which is centered near x — 1/2, the "CZ wave function" has a "double humped" form, with maxima near 



the extremes of £. It has been the subject of much controversy, as will be discussed in Sec. |3.1 
0.30 
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Figure 4: CZ wave function (solid) and asymptotic wave function (dashed). 

P 2| , we first perform a formal Taylor expansion 



To derive sum rules for moments of the wave function 
of the quark field u(z) in Eq. (p2|), 



•) = E 



(d + ) n u(0) . 



(25) 



n=0 



Substituting the resulting expression into Eq. (|23|), and carrying out the z integrals, we derive the following 
expansion in local operators, 



E 



n=0 " ! P 



, ( ^ (°l^ + 75 (*r«(o)k(p)> 



Moments of W with respect to x then pick out individual matrix elements: 



(p + ) n+1 / dxx"^(x,/i 2 ) = i n (0|J„(0)|7r(p)> , 
Jo 



where 



J„(0) = J(0) 7 + 7 5(5 + ) n u(0). 



(26) 



(27) 



(28) 



Analogous relations between moments of a light-cone distribution and matrix elements of local operators are 
familiar from DIS. 

We now consider the specific Green function (correlator) 

G n {p 2 ,p + ) = I d 4 zc^(0|T[J n (z)J (0)] |0) , 



(29) 

with J n {x) given by Eq. (|28|). Such a two-field Green function enjoys the analyticity structure shown in Fig. ||. 
At fixed p + > 0, G n (p 2 ,p^) is an analytic function in the complex p 2 plane, except for poles and branch cuts 
along the real, positive p 2 axis. By Cauchy's theorem, the integral of G n (p 2 ,p + ) along the two contours Ca 
and Cb in Fig. |^ give the same result. The value p 2 = sq where these contours meet is sometimes called the 
"duality interval" , referring to the complementary (dual) manners in which they are evaluated. 

Contour Ca is evaluated using our knowledge of hadron spectroscopy. Because the contour runs around 
the real axis, the integral is given by the imaginary part of G n (p 2 ,p + ) - a sum of delta function contributions 
from hadronic bound states with the quantum numbers of the pion, plus possible multiparticle continuum 
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contributions. To emphasize the lowest-lying state, in this case the pion, we multiply G n (p 2 ,p + ) by the "entire" 
function e~ p ' M , with M 2 an adjustable mass. Then, a short calculation for the integral along Ca gives 

I A {M 2 ) = I J^ Gn (p 2 ,P + )e-" 2 / M2 
Jc A 2ttiM^ 

(i P +r +i 



M 



r2 



x n ){x°)+R A , (30) 



where (x n ) is the weighted integral of in Eq. (|27|), and Ra is a remainder, associated with higher-mass 
resonances such as the A\, and the continuum. This integral is referred to as a Borel transform. 

The corresponding integral Ib{M 2 ) taken along contour Cb is computed quite differently. Along Cb, the 
integrand is evaluated far from any resonances, and we may hope that it behaves as it does for Euclidean 
p 2 < 0, where the OPE applies. Rules for its calculation are straightforward but rather technical. The basic 
structure of the answer, however, is readily expressed as a sum of a perturbative term plus two nonperturbative 
contributions, from the gluon condensate {G flu {A)G^ u (A)) and the quark condensate (qq)o, 

I B (M 2 ) = h { I n) (M 2 ,p+) + h^(G 2 ) + h\f q)2 ((qq) a ) 2 + R B , (31) 

where Rb represents corrections. All the h's ("coefficient function of the OPE) are computed in perturbation 
theory. 

The values of M 2 and so are to be chosen to minimize Ra in Eq. (|30"1). Values of < qq >o and < G 2 >o 
may be found from the analysis of e + e~ — > hadrons [p9| . Finally the coefficient functions depend on a 
renormalization scale /jq 2 . Combining these choices and parameters, and setting I a = Ib, we may therefore 
determine (x n ), or equivalently (£"), with £ = 1 — 2x the relative fractional momentum. The CZ wave function 



in Eq. (24) above was found by fitting its moments to those found by the sum rules. 



2.4 Beyond the Pion 

2.4.1 Generalizations Most of the developments outlined above for the pion apply as well to electromagnetic 
form factors for other hadrons, especially baryons |3^, ^|, |3^] and also resonance production, as well as vector 
mesons and kaons || [32[ [}7| [58|. The form factors of baryons are determined by three-quark valence wave 
functions, and for both vector mesons and baryons nontrivial spin structure must be taken into account. So 
long as transverse degrees of freedom may be neglected, however, spin may be described in terms of conserved 
hclicity, where the helicity of a hadron is given by the sum of the helicities of the its partons. A PQCD treatment 
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of violations of helicity conservation has been proposed in [[L6| |17| . We shall have occasion below to review some 
of the successes and limitations of this rich constellation of predictions for hadronic form factors. 

For example, the wave function of a proton is a sum of terms describing total helicity ±1/2, times functions 
<pi(xi, X2, %3, fJ. 2 ) with J2i x i = The application of evolution analysis to these wave functions shows that 
asymptotically that have the simple form, 

<pASY( x i) — const, x X1X2X3 . (32) 

In this case, no readily observed decay amplitude is available to normalize the asymptotic wave function, and 
hence the proton's form factor. A Sudakov analysis of the proton wave function and form factor is also possible, 
with the same general properties as for the pion. [Q, It involves two transverse separations, however, and is 
correspondingly more complex. 

Another important difference between the proton and the pion is in the baryonic analogue of Eq. ([?]) for 
helicity form factors G (see below), which we may represent schematically as 

G(Q 2 ) — J dxidx 2 dyidi/2 (/} P (yj,fi 2 )T G (yj,Xi,Q) (/) P (xi,fi 2 ) , (33) 

where Tq, and hence G, is proportional to Q~ 4 , and begins at order a 2 . Here, in contrast to Eq. (Q) for 
the pion, however, the perturbative expansion of the "hard-scattering" function Tq receives infrared divergent 
contributions from regions that resemble the Feynman mechanism, in which one quark carries essentially all 
of the proton's momentum, beginning at two loop corrections j^, [l(|. Such regions are suppressed by Sudakov 
corrections. Progress has been made in quantifying this observation for valence PQCD, by introducing transverse 
degrees of freedom for baryons, as for the pion, but a complete formalism for baryon form factors, even to leading 
power in Q 2 , remains for the future. 

2.4.2 Baryon Helicity Matrix Elements. For use below, let us define electromagnetic helicity matrix 
elements for nucleons. Taking into account resonance production, an initial state with helicity A = 1/2 may 
become a final state with A' = 1/2 or 3/2. Transitions between a nucleon state \N >, and final state \N' > can 
be expressed in terms of dimensionless helicity matrix elements, 



G H = — - <N',\'\e»- J M \N,l/2>. (34) 

ZM jv 

This notation follows The polarization vectors e ±,Q correspond to right (G+) and left (G_) circularly 

polarized photons and longitudinally (Go) polarized photons, respectively. G+,Go,G_ describe transitions 
in which AA = 0, 1 and 2, respectively. Assuming that helicity is conserved, valence PQCD suggests that 
G+ cx QGo oc Q 2 G-. For elastic scattering, since the recoil nucleon has spin 1/2, only the helicity conserving 
G+ (AA = ) and non-conserving Go (AA = 0) contribute. 



2.5 Nonasymptotic Form Factors 

The valence PQCD results above determine the form factor at very high Q 2 . How high one must be, however, is 
a matter of debate (see below). It is therefore important to develop treatments of the transition to asymptotic 
behavior. The evolution of wave functions is a step in this direction, but at moderate Q 2 , it is necessary to 
apply methods, or develop models that take into account processes that are suppressed even by powers of Q 2 
at high energy. These include the soft processes discussed above. 

2.5.1 Sum rules for Form Factors. Refs. [[Io| and pjj ], have utilized the sum rule approach to directly 
obtain form factors, without the intermediate step of determining wave functions. This approach, as described 
above, depends on the analyticity properties of Green functions that are associated with form factors. It is 
thus not a dynamical theory of soft or hard interactions, but relies on general properties of QCD, such as the 
OPE, in addition to perturbative calculations. For a hybrid approach, with features of both QCD sum rules 
and valence PQCD, see Q. 

For the pion form factor the relevant Green function may be expressed as 



T^ Pl ,p 2 ) = i 2 J e^-y + ^ z {0\T[J{y)Jl m {0)J{z)\Q)d 4 zd i y, 



(35) 
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with J = Jq defined as in Eq. (|2£]), and J em the electromagnetic current. In terms of a related scalar amplitude 
T, this Green function possesses a double dispersion relation, 



1 pQO pQO { 

2 2^ 1 / j_ / j„ P(Sl,S2, 



T(p-L,p 2 ,Q ) = Zz I <M ds 2 {si _ pl) ^ s2 _ ply (36) 



The spectral function contains a pion pole which defines the pion form factor, p 7nr (si,S2, Q ) = 2-n - 2 f 2 F^{Q 2 )5{s\- 
m 2 r )S(s2 — m 2 ), as well as a continuum above the 3-pion threshold, which also includes the broad A\ state. The 
form factor F n is extracted by relating the two contours of Fig. 0, this time in both variables si and S2- In 
the Borel transform is replaced by a simple integral (Af 2 — > oo), and so ~ 8-7T 2 / 2 ~ 0.7GeV 2 is adjusted 
to reflect this choice, known as "local duality". This leads to a relation between F 7T (Q 2 ) and the lowest-order 
perturbative contribution to p, which may be evaluated to give 

F.m-i- (1+6s ° /g2 > ■ (37) 

(l + 4s„/<3 2 ) 3 ' 

Expanding in inverse powers of Q, this expression behaves as Q for large momentum transfers, and is thus 
eventually nonleading compared to the perturbative prediction (|2(]). Nevertheless, as we shall see, it gives a 
viable fit to the available data, which implies at the least that "soft physics" plays an important role in the 
charged pion form factor at present energies. Beyond lowest order, p includes gluonic corrections, which appear 
to correspond to the hard gluons of valence QCD. Similar methods have be used to treat baryon form factors 

ilil- 

2.5.2 Models. Unfortunately the complexity of soft processes in QCD does not lend them to simple physical 
models. Their description in terms of fundamental QCD is one of the outstanding theoretical challenges in 
the theory. There have, however, been useful attempts to bridge the low and high Q 2 regions with various 
phcnomcnological or empirical approaches, concentrating on nucleon form factors. 

The generalized vector dominance model (VDM) or hybrid model of |i"5| begins with the VDM, which yields 
the requisite low Q 2 form-factor. Additional terms join VDM form-factors smoothly to PQCD expectations at 
high Q 2 (Gm oc Q~ 4 and Ge oc Q~ 6 ). With the appropriate choice of parameters an excellent agreement with 
the G£f data is achieved over the entire range of available Q 2 . Agreement with the other elastic form factors, 
however, turns out to be poor in the light of more recent data. 

The Constituent quark model has been been modified, and relativized to extend their validity into the few 
GeV 2 region of Q 2 [ll| . For example, in the calculation of hadronic form factors in jl9| , the constituent 

quarks, of mass ~ .33 GeV, have wave functions which are solutions to a potential derived from a quark-quark 
interaction model. In a light cone frame the wave function takes the form i/j(x,pt) ~ X(x)P(x,pt). The range 
of pt in the model wave function effectively has an ultraviolet cutoff so that the one-gluon perturbative parts are 
not included in the derived form factor. With reasonable choice of X the soft components play an important, 
and even dominant role over the entire range of measured Q 2 . However, there models are not rigorous enough 
to make precise predictions. 

The diquark model [ [i"8[ [i"9| , |50|| assumes that the baryon distribution function can be expressed in terms of 
two constituents, a quark and a diquark, which consists of a correlated quark pair. The diquark structure allows 
for helicity non-conservation, and thus at some level can also account for soft processes. The diquark becomes 
completely equivalent to the valence PQCD model in the high Q 2 limit. Its several parameters can be tuned to 
give a good fit over the entire range of G£f, including the transition Q 2 range. 



3 Experimental Status of Hadronic Form Factors 
3.1 Pion Form Factors 

In this section we will discuss the ir + and ir° form factors as obtained in the reactions p(e, e'it + )n and e + + e~ — * 
7T°, respectively. Given the relative simplicty of the mesonic valence state, we might expect perturbative analysis 
to apply at lower momentum transfers for pions than for nucleons. We discuss the successes and shortcomings 
of the valence QCD approach in explaining the data, and also point out important uncertainties in the data 
itself at high Q 2 . 

3.1.1 The Charged Pion Form Factor. The tt + form-factor is obtained by studying electroproduction 
on a hydrogen target (see Fig. ||). The aim is to separate the "i-channel process", in which the electron scatters 
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from a nearly on-shcll virtual pion emitted from the proton. This t-channel cross section, which is due to the 
exchange of a longitudinal (L) photon, determines the pion form factor, though the relation 



f 2 AQ 2 ), 



(t — m 2 ) 2 

where t is the squared momentum transfer to the nucleon, and <? 2 jVArW * s ^ ne ^NN coupling. 



(38) 





SnNh 
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Figure 6: a) Interpretation of t-channcl n + production in terms of the pion form factor F % . b) Lowest order 
diagram for the tt° form factor. 



Nearly all the existing high Q 2 data, shown in Fig. [3, were obtained at Cornell |5l], 52 , ^3|. Care, however, 
must be exercized in the interpretation of the higher Q 1 points, which do not include systematic errors. The 
reason for this uncertainty is that the separation of (7l from the complete cross section requires measurements 
at different electron scattering angles at the same Q 2 . This "Rosenbluth separation", was not practical at the 
highest Q 2 in this experiment. For Q 2 > 4 GeV 2 the (unwanted) transverse cross section was estimated from 
an extrapolation of low Q 2 data, and subtracted by hand. Thus, although reliable data exist for Q 2 < 3 GeV 2 , 
the 6.3 and 9.7 GeV 2 points provide little help in distinguishing between theoretical models. 




Figure 7: The form factor of the 7r + meson Q 2 F^{Q 2 ) vs Q 2 . The data are from |H| ^2], |53|. See the text for 
comments on the interpretation of the higher Q 2 data. The dot-dash curves labelled CZ and ASY are obtained 
from Eq. (39), using (j)cz and 4>asy respectively. The dashed curve labelled SR is the direct sum rule result 
of The solid curve LS is from ||, using the CZ valence quark distribution amplitude, and including the 
effects of Sudakov suppression. 

There are also important theoretical issues in the extraction of the data. For instance, the struck pion is 
off-shell, and one must extrapolate to the physical pion pole at t — —m 2 . Uncertainties in the t dependence 
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of g 2 NN (t) also lead to uncertainties in F n . In addition, the reliability of high Q 2 form factors extracted in 
this manner has been questioned by |54|] , who claim that other hard, non-resonant processes compete with the 
i-channel process, and may be difficult to separate from it. These objections aside, an important future goal is 
to extend the pion form factor data to higher Q 2 f35| . 

3.1.2 Comparison with Theory. In the valence PQCD framework, the pion form factor may be written 
in factorized form as in Eq. (f7|). Treating the hard-scattering at lowest order, with f n ~ 93 MeV, we have 

FAQ 2 )== l6nCl Z s{K2) \I\ 2 with J=/ 1 ^^M. (39) 



Q 2 

This formula, with a valence quark distribution amplitude derived from QCD sum rules fl32|| , denoted 4>cz, 
gives a pion form factor in rough agreement with the data as shown in Fig. [?]. In obtaining this, the variation 
in as was fit to the evaluated data Q, with n 2 = Q 2 /4. The asymptotic distribution amplitude 4>asy, Eq. 
(|l9|), which yields Eq. (20), seriously underestimates the data. Refering to Fig. |[ the difference is that 4>cz, 
Eq. (HJ), has a "double- hump" structure, concentrated near x ~ and 1, and hence yields a larger value for 
/ than the more central 4>asy- This apparent success inspired many theoretical papers based upon the "QCD 
sum rule" technique for describing exclusive reactions. The authors of Jl^, on the other hand, observed that 
with <j)cz, Eq. (|9]) is dominated by soft gluon momenta k 2 (— xyQ 2 ), near the end-point regions discussed 
above. They argued that Eq. (|3^), or for that matter PQCD, is invalid in the kinematic regime where data is 
available, because higher-order perturbative corrections would be uncontrollably large for gluons of such low 
momenta. If one cuts off the integral in Eq. ( |39| ) below a minimum gluon invariant mass, say k 2 ~ .5 GeV 2 , 
one derives a much smaller "legal" part of the form factor ( ~ 10 - 20 percent remains for Q 2 between 5 and 10 
GeV 2 ). 

Roughly, proponents of valence PQCD were faced with the dual problems of how to keep the main contri- 
butions to the x integral in Eq. ( |39"| ) away from the cndpoints, at the same time enhancing their values relative 
to the simple use of (j^ASY- One way of doing this is to resum a selection of higher-order corrections into the 
argument of the strong coupling. Choosing [i 2 — xyQ 2 in T in (^) results in a significant enhancement, because 
the perturbative running coupling grows as its scale decreases. This running coupling, however, diverges for 
x yQ 2 — Aq CD , which requires the introduction of a scale below which the coupling is "frozen" . The result 
is naturally quite sensitive to the cutoff, but it can give a reasonable result without dipping too far into the 
nonperturbative region |57j . 

In a related development, it was argued that transverse degrees of freedom should not be neglected, and 
indeed mimic a gluon effective mass, which suppresses the blowup near x, y = |Q. We have already seen 
how Sudakov resummation of transverse degrees of freedom in Eq. (|2l]) results in a naturally self-consistent 
calculation of the form factor, without cutoffs ||. In this case, the CZ distribution amplitudes continued to 
account for the existing data, when the enhancement associated with the running coupling was included. These 
results for the n + form factor are plotted in Fig. \q. 

The calculation of Q has been generalized in [[59[ , who specifically included an "intrinsic" transverse wave 
function. That is, in Eq. fl2l| ) above, they replaced exp(— S) — > exp(-S') T,(x, b). Using a model, Gaussian shape 
for S, they found that this further protected the resulting form factor from the soft region, but also further 
supressed the hard part of the form factor below the data. Of couse, this proceedure introduced an additional 
parameter, in the Gaussian, and it included a constituent quark mass. 



In an alternatative approach (described in Sec. |2.5[ l above) the direct prediction of F n from QCD sum rules, 
Eq. (|7]) 0, accounts for most of the measured F n) without including gluon exchange into its perturbative 
calculation, even though the resulting expression decays as Q 4 a t higher Q 2 . In this and other alternatives 
to the valence quark picture, the apparent scaling with Q 2 of the present data is interpreted as something 
of an accident. The extra contribution of a gluon exchanged between quarks, which produces Q~ 2 behavior 
asymptotically, has been estimated and leads to a modest increase at the highest available Q 2 . These 
two contributions are referred to as "soft" and "hard" ]6(J], the latter being identified with the valence PQCD 
prediction. 

Other publications continue to focus on the relative importance of soft and hard processes, and in particular 
how to deal with the difficult soft sector. Examples are |^2|, |3l], [H], |5^, |3Cj, who all conclude that soft 
processes are important for Q 2 corresponding to the existing data. 

Timelike (s = q 2 > 0) Pion Form Factor. This is obtained in the reaction e + + e~ -^7* — > tt + + tc~ . 
Only one data point exists in the multi-GeV 2 region, at s — Mju ~ 9.6 GeV 2 , obtained from the ratio 

(J/ip — * 7r + 7r~)/(J/?/> — * e + e~) by |32|. This point appears to be more reliable than the higher Q 2 space-like 
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data. Its calculation is identical to the space- like case in most respects, and it would be useful to obtain timelike 
data over a range of Q 2 . This process was calculated Q with valence PQCD techniques, employing evolution 
and Sudakov suppression, as in Q. The ratio of experimental timelike to spacelike form factors, about 2, is 
consistent with the valence PQCD calculation, although the overall normalization is low by a factor of two or 
more, depending on the light-cone distribution amplitude employed. 

3.1.3 The 7T° Form Factor. The 7 + 7* — > 7r° form factor is expected to be a particularly good test for the 
pion's valence distribution amplitude, since at lowest order in the hard scattering it is a pure QED processes 
(see Fig. ^b). Higher Fock state contributions are suppressed by powers of a s (Q 2 )/Q 2 . In addition, there is no 
analogue of the "soft", Feynman mechanism contributions, which require an incoming and an outgoing pion. 

Experimentally, the 7r° form factor can be studied via either the Primakoff effect or virtual Compton scatter- 
ing. The former is accessible in e + — e~ colliders, while the latter is more appropriate to fixed target machines. 
Fig. H includes data of the CLEO-II group, J64| which reported measurements up to Q 2 ^ 8 GeV 2 , from reactions 
e + e" — ► 7r + X. 



0.30 




Figure 8: Form factor in the few GeV 2 range for 7 + 7* — * 7r°. Data for Q 2 > 3 GeV 2 are from |j4j, and 
lower Q 2 data from |6(|. The dot-dash curve labelled CZ is the result of using tficz in the integral / in Eq. 
(40), and the solid curve labelled ASY is obtained when using (j>ASY- To obtain these curves, Eq. (40) has been 
modified so that the form factor joins smoothly with the known value at Q 2 — 0, using a generalization of the 
prescription of M. The horizontal line at the right labelled ASY-limit is the high Q 2 limit, using 4>asy- The 
dashed curve labelled SR is the soft form factor obtained directly from QCD sum rules |57J. 



Working to lowest (zeroth) order in a s , by analogy with Eq. (p9|), the relationship between the ir° form 
factor and the valence quark distribution amplitudes is 

F 77 .o(Q 2 ) = -^L_/, I = J\ X *M, (40) 

A calculation [65| following Q, including Sudakov effects, and an intrinsic transverse distribution amplitude as 
above, 

dx— ^(xjEfo b)f H (x, b, Q)e' s ^- b ^ , (41) 

appears to account well for the tt° form factor. The results with 0^ = 4>cz and 4>asy are shown in Fig. ||. To 
compare theory with experiment in this lower Q 2 region, Eq. ( (40|) has been modified so that the form factor joins 
smoothly with the known value at Q 2 — 0, using a generalization of the prescription of [^3| . In this case, <\>asy 
accounts for the data, while <&cz overshoots it. A recent sum rule computation of this form factor (6^] also 
accounts well for the data. Because the hard-scattering in the 7r° form factor starts, as in Eq. (fh)]), at zeroth 
order in a s , the sum rule and valence PQCD approaches both begin with the same perturbative diagrams. 
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Finally, one may determine the integral / directly by fitting Eq. ( f4(i| ) to the ir° data, as shown in Fig. |[ The 
result is close to the value of / for 4>asy- 

Comparing the it and tt + form factors, the success of the asymptotic distribution amplitude in the former 
suggests that 4>asy should be used to compute the the valence PQCD contribution to the latter as well. Then, 
however, the valence contribution with lowest order gluon exchange accounts for less than one half of the ir + data 
(see, for instance, Q). We conclude that, if the charged pion data is accurate at all, either non- valence (soft) 
contributions, or higher-order contributions in valence PQCD, must play an important role. We note a lack of 
need for the soft mechanism in the ir° form factor, which is consistent with these observations. Higher-order 
hard corrections are also different in the two form factors, however, so it is difficult to draw a final conclusion 
without further study. We consider, however, that it is likely that the soft mechanism plays an important, and 
possibly dominant, role in the region of a few GeV 2 where reliable data exist. 



3.2 Nucleon Form Factors 

In this section and the next we consider nucleon elastic form factors and transition form factors involving 
resonant states of nucleons. Since there exist stable on-shell baryon targets, the nucleons, and there are a large 
variety of final states of spin and isospin, the resonances, a wealth of experimental information can be accessed. 
Nevertheless, rather limited data exist at high Q 2 . 

3.2.1 Nucleon Elastic Form Factors. The elastic electron-nucleon cross section, expressed in terms of 
the Sachs form factors, is 



da 



= cTMfrec ' - 1 + T + 2r|G M | 2 tan* 0/2 j . (42) 

Here, cm is the Mott cross section for scattering from a point object, f rec = E'/E is a recoil factor, r = 
Q 2 /4Mfj and k is the anomalous magnetic moment of the nucleon, in nuclear magnetons: up = 1.79, and 
k n = -1.91. 



The helicity matrix elements G±,o, defined above in Sec. f^J, are related to the Sachs form factors, and the 
Fermi (Pi) and Pauli form factors (F2) as follows: 

Go = Ge = *i - jjjj-itfi. («) 

For elastic scattering from nucleons there are two helicity conserving and two helicity non-conserving form 
factors, G£f, G^ and G^, G^, respectively. At low Q 2 (less than one or two GeV 2 ) all the form factors are 
consistent with a dipole Q 2 dependence, 1/(1 + Q 2 /M 2 ) 2 , with M ~ 0.71 GeV 2 . 

At high Q 2 , valence PQCD predicts that the Q 2 behavior of the helicity conserving form factors G^ or 
GjJJ- should follow a 2 (Q 2 )/Q 4 (see Sec. [2.4| .l), where we recall that a s {Q 2 ) decreases logarithmically in Q 2 . In 
addition, their magnitudes are determined by relations like ( |33| ) using nucleon wave functions, either of the 
asymptotic form (Eq. ([52])), or as found, for instance, from sum rules. The helicity non-conserving form factors 
should Gv N , should fall as G^f /Q 2 §. In the high Q 2 limit of Eq. @, Pi ~ G M oc G+. 

Fig. fl summarizes what is known experimentally about these four form factors, which were mostly obtained 
at SLAC. G£f is known best, followed in order by G^, G^ and G^. We will consider each in turn. 

3.2.2 The Proton Magnetic Form Factor. Only G£f has been measured at high Q 2 . In the low Q 2 limit 
G^ and G^ are comparable, and can be separated with comparable accuracy by a "Rosenbluth" separation. 
Separated form factors only exist out to Q 2 ~ 9 GeV 2 and unseparated data exist up to Q 2 ~ 31 GeV 2 p9[ |. 
However, at lower Q 2 it is observed that G^ is much smaller than G^, and that they are roughly proportional. 
Since at higher Q 2 the G^ contribution is kinematically suppressed (see Eq. Ref. (fl^)), p9| | estimated G^, 
assuming only that G^ does not grow anomalously. The result is presented as a measurement of G£f. 

Much theoretical work has focussed on the application of the valence PQCD techniques described above to 
the calculation of the helicity conserving G£f (actually Pi) |52|, 57 . The broad issues are similar to 



those discussed above for the pion. An advantage relative to the pion, however, is that, because we can scatter 
electrons from on-shell protons, G£f is relatively unambiguous over a larger range of Q 2 . A disadvantage is that 
the proton's three valence quarks make it theoretically more complex. 

It was observed quite early that the Q 2 dependence of G^ is in agreement with quark counting (and hence 
valence PQCD) predictions. We have already encountered the basic methods and arguments in our discussion 
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Figure 9: The clastic form factors as a function of Q 2 divided by the dipole shape: a)G^ (spacelike): • - pCf , 
- [pi |l|. b)Gg. The curves are: solid f|, dashed ||, dot-dash Jp. G P M (timelike): • - - Q. h)G p E : 
• |70| . The curves are: dashed ^o|, dot-dash jl5|]. c)Q 2 F2/Fi: • |7(| . The curves are: dashed |2Q] , dot-dash - 
|5J. d)G^: • ||, l7|, x (7l. The curves are: solid $7\, dashed M, dot-dash @. e)Gf : • x @. 
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for the charged pion form factor. Once again, calculations based on lowest-order gluon exchange and asymptotic 
distribution amplitudes fall far below the data. After the initial development of sum rule distribution amplitudes 
[ [32| , however, the situation appeared to improve. This was taken as compelling evidence of the applicability of 
valence PQCD techniques at measurable Q 2 . This conclusion has been the focus of many papers, sometimes 
quite contentious, by both proponents and detractors of the use of valence PQCD at accessible Q 2 . 

An example is the calculation of J57|, based on wave functions derived from QCD sum rules (Sec. 
The result of using (f>Q Z from J32| is plotted along with the data in Fig. ^a. Once again, the coupling is forced 
to run with the virtuality of the exchanged gluons, down to a mass scale at which it is frozen. Good fits to the 
data were obtained, when this scale (termed an effective gluon mass) is 0.3 GeV. The curve in Fig. ^a is not 
applicable at low Q 2 , since it is based on leading order PQCD. 

This approach has been strongly criticized [jl9|, ^0), for the proton as for the pion. The basic question is 
whether the major contribution to the form factor comes from gluon exchange at low virtuality, where higher- 
order contributions are not under control. Indeed, sum rules seem to suggest asymmetric nucleon distribution 
amplitudes, which, as in the case of the pion, enhance contributions from low gluon virtuality. At the very least, 
this produces strong sensitivity to the mass at which the coupling is frozen, and shakes our confidence in the 
self-consistency of the valence calculation. Also, it was concluded |7(| f7?J that the uncertainties in obtaining 
reliable distribution functions from sum rules, given the experimental uncertainties in the condensates, are so 
great that the distribution amplitudes are essentially undetermined from sum rules alone. We may also note 
that the lattice calculation of [ p6| supports a rather symmetric wave function in the nucleon. 

Finally, as for the charged pion form factor, the inclusion of transverse momentum effects |5^, pi stabilizes 
valence PQCD calculations and improves their self-consistency, while generally reducing them. Thus, Ref. [78], 
following up on the pion calculation of ]^|, recalculated F\ (or Gj^-) according to the techniques of jq, 36], 
including intrinsic k± componants in the distribution amplitude, (f>p(x) — ► <pp(x)Qp(x,b), with Q a Gaussian. 
In contrast to the pion form factor, the Sudakov resummation for the proton form factor leaves a sensitivity 
to large b in a corner of the b, x space, necessitating the inclusion of an "infrared" cutoff, that is, a maximum 
transverse separation in the distribution amplitude. It should be noted, however, that alternate resummations 
that suppress all large b should be possible, although they have not been explored in the literature. As in the pion 
case, the extended calculation of |79] reduces the hard scattering form factor significantly below experimental 
data, for both 4>q Z and 4>asy- 

G£f at time-like momentum transfer can extend the range of Q 2 and, together with space-like data, can 
further constrain theory. The timelike proton form factor has been measured for three values of Q 2 near 10 
GeV 2 [f79|. This, with lower Q 2 data j8(| are also shown in Fig. ||a. As in the pion case, a factor of about two in 
the ratio for the space-like and time-like form factors is consistent with expectations from valence PQCD p3| . 

3.2.3 Proton Electric Form Factor. As Q 2 increases, kinematic suppression of the contribution of in 
Eq. ([l2]) makes a Rosenbluth separation less and less accurate. As a result, once Q 2 ~ few GeV 2 , the errors on 
available data for G|J are significantly worse than for G£f. The most recent data |nj obtained by Rosenbluth 
separation exhibits much smaller errors than previous data, and extends the measured range out to Q ~ 9 
GeV 2 . The data, shown in Fig. follow a dipole shape over the entire range of Q 2 to within the limited 
accuracy. 

As indicated above, because G^ oc Gq , which is helicity non conserving, at high Q 2 the ratio Q 2 G^/G^j 
or Q 2 F^ /Ff should approach a constant. As seen in Fig. ^c, it appears to do so. This qualitative success of 
valence PQCD in the 5-10 GeV 2 range makes it attractive to extend the experimental range of accurate G^ 
data to higher values of Q 2 . A decrease in the ratio Q 2 G 1 ^/G 1 ^ I for large Q 2 of 20 GeV 2 , say, might be a signal 
that soft processes are still dominant over hard processes in this range. 

Of course, at increasing Q 2 the Rosenbluth separation becomes more difficult. Other methods, involving 
polarized beam and target or polarized beam and proton recoil polarimeter ||S2| , |83| , [S4|| , which measure the ratio 
G^/G^j, become more favorable. Using such techniques, it will be possible to extend measurements of G|J to 
higher Q 2 . 

3.2.4 Neutron Form Factors. Form factors of neutrons are difficult to obtain, because there are no free 
neutron targets. Most of the available data were obtained in quasielastic scattering from deuterons, in which 
the proton contribution is subtracted. This method has intrinsic uncertainies, since one must deconvolute from 
the quasielastic peak the contributing neutron and proton nuclear wave functions, which must be independently 
known, as well as the intrusive tails of the inelastic processes, which are also broadened by Fermi motion. 
This becomes increasingly difficult with increasing Q 2 , as the contribution from quasielastic scattering relative 
to the inelastic processes decreases dramatically. Eventually, the tail of the inelastic background dominates, 
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and the the extraction of the quasielastic peak becomes extremely sensitive to uncertainties in the modeling of 
inelastic processes. Thus, at this time data on G^ and G^ are limited to the range Q 2 < 10 GcV 2 and 4 GeV 2 
respectively. 

There are various ways of improving the situation. The detection in coincidence of the struck neutron 
along with the electron can effectively eliminate the quasielastic proton contribution, and significantly reduce 
background due to inelastic processes. For G^, which is much smaller than G^, polarization asymmetry 
techniques can yield the ratio G^/G^i- This method has been employed successfully at lower Q 2 pq ], and is 
currently planned |36, 87 for the few GeV 2 range. Neutron form factor data in the GeV 2 region were obtained 
at SLAC J72j, emploing careful Rosenbluth L/T separations of single arm cross section measurements. The 
available data are summarized in Fig. ^c and Fig. ^d. The G^ data are consistent with the dipole shape over 
the entire range of Q 2 , athough there is significant variance between data sets below about 1 GeV 2 . The data 
for G^ is consistent with zero up to the highest Q 2 , although the errors are quite large. 

The value of this data, even though the range is mostly limited to the region where soft processes may still 
dominate, is quite apparent. All the theoretical non-valence PQCD curves deviate from the data on G^ with 
increasing Q 2 . Examples shown are the hybrid |45| and the QCD sum rule result of (20). The constituent quark 
|p7| , and vector dominance |pq| models also appear to diverge monotonically with increasing Q 2 . The data on 
G^ clearly eliminates the hybrid model, whereas the VDM and QCD sum rule based calculations are consistant 
with zero over the Q 2 range. 

To make further use of the selectivity of the nucleon form factors, it will be important to obtain data on 
Gj, G^ and G 1 ^. at Q 2 greater than the present limits. Such experiments for G^ and G^ have been proposed 



for future facilities |89, Q. For G^ at least, its scaling as Q~ 4 over such a large range suggests that valence 
PQCD is relevant to its description. The soft mechanism may, however, also play an important role, especially 
at moderate Q 2 . The clarification of this role is an important project for theory and experiment. 

Only global tests involving all available form factors can hope to seriously select among varying points of 
view. We stress the importance of measuring the helicity non-conserving form factors to as high Q 2 as possible, 
since in valence PQCD they are driven by non-leading processes, and therefore offer important constraints on 
the relative importance of soft and hard processes with varying Q 2 . 



3.3 Baryon Resonance Amplitudes and Form Factors 

The study of transition form-factors to excited baryons at high Q 2 can make an important contribution to our 
knowledge of hadronic structure. Fig. [To] shows the virtual photon cross section at Q 2 = 1 GeV 2 as a function 
of baryon invariant mass W. For W < 2 GeV, the most significant feature is the existence of three maxima, 
known as the first, second and third resonance regions. In this interval there are about 20 known resonances. 
These are denoted -£27,2,7 (WO, where L is the angular momentum of the single pion decay, and I and J are 
respectively the resonance isospin and spin. However, except for the first, which is due to the A(1232) the 
resonances are largely overlapping, even with a significant non-resonant underlay. In future programs, the 
separation of the contributing electromagnetic multipoles will require measurement of exclusive reactions such 
as (e, eV) and (e, e'77) to as high Q 2 as possible, with polarized beams and targets. 

The second resonance region is dominated by two strong negative parity states, the Di3(1520) and the 
Sn(1535). At low Q 2 ( < 1 GeV 2 ) the D 13 (1520) is dominant, whereas at higher Q 2 ( > 3 GeV 2 ) the 5u(1535) 
dominates. The Roper resonance, the Pn(1440), has not yet been definitely observed at Q 2 > 0, but is of 
considerable interest since there is speculation regarding its character |^|. In the third resonance region, the 
largest excitation at low Q 2 is the Fis(1680). The relative strength of the other states is not well determined, 
especially at increasing Q 2 . At low Q 2 the excitations indicated in Fig. |lO|have been rather successfully described 
in terms of the constituent quark model. 

The current experimental situation is that exclusive (e, e'lr) and (e, e' rf) data exist only up to Q 2 = 4 
GeV 2 . Although there is a total absence of exclusive data above Q 2 =3 GeV 2 , there are inclusive data in the 
resonance region obtained mostly at SLAC (see references in |92]j)- Although the statistical accuracy becomes 
poor at high Q 2 , the three peaks near W — 1232, 1535 and 1680 MeV remain prominent, with the A(1232) 
obviously decreasing with increasing Q 2 relative to the other two. After subtraction of phenomenological non- 
resonant backgrounds the peaks were fit with resonance functions ( [p2| ) to extract transverse form factors, 
|G T (Q 2 )| 2 = (|G+| 2 + |G_| 2 )/2r, where r = Q 2 /4M 2 . 

The form factors, are shown in Fig. [□], relative to a dipole shape. Also shown at lower Q 2 are form factors 
extracted from data obtained earlier from exclusive (e,e',p)ir° and (e, e',p)r] experiments. 
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Figure 10: The virtual photon cross section for proton excitation at Q 2 = 1 GeV 2 . The data are reconstructed 
from an evaluation by |p0[ . The known contributing states are indicated at the bottom by vertical lines. The 
largest contributing states are indicated by long vertical lines. 



Fig. |Tl| shows that the form factors obtained for the second and third resonance regions are consistent with 
a Q~ 4 dependence, although with large statistical uncertainty. On the other hand, the A(1232) form factor is 
decreasing relative to both the elastic as well as the second and third resonances. Since this result is obtained 
from inclusive data there are systematic uncertainties in the extraction |Q , and more recent analysis of the 
available data pjj] indicates the extent in Q 2 of the decrease is yet resolved. 

3.3.2 A(1232): Transition Multipoles Since the A(1232) has J = 3/2 there are three contributing 
multipoles, Ei+,Mi+, and S\+ whose relative contributions are model dependent. Thus, this is a favorable 
case for studying models of baryon structure. At low Q 2 in a pure SU (6) non-relativistic CQM the N — > A 
transition is purely Af 1+ in character, involving a single-quark spin-flip with AL = 0. An E\ + contribution 
is not permitted, since the A and N are both in L = states, which cannot be connected by an operator 
involving L > 0. The addition of a residual quark-quark color magnetic interaction adds higher L components 
to the A wave function, and thus introduces a small E\+ component, of perhaps a few percent. At Q 2 — the 
experimental data supports the constituent quark model prediction of Mi+ dominance extremely well. Recent 
data Q bears this out. The data from from Mainz @ reports a ratio E 1+ /M 1+ = - .025 ± .002 ± .002. This 
ratio remains very small up to Q 2 about 1 GeV 2 , beyond which there is very little data. There exist some 
earlier data at Q 2 = 3 GeV 2 [p7| , which has been evaluated by Q, suggesting that E\ + /Mi + is increasing, but 
with large errors, Re(Ei + /Mi + ) — 0.06 ± 0.02 ± .03, and we must conclude that the magnitude of E\ + /M\ + at 
Q 2 — 3 GeV 2 remains uncertain. Recently (9i| exclusive data were obtained at CEBAF at Q 2 = 3 and 4 GeV 2 , 
for the A(1232) and 5n(1535), but at the time of writing the analysis is not complete. 

At high Q 2 , valence PQCD predicts that only helicity concerving amplitudes should contribute. The multi- 
pole amplitudes for single pion production may be expressed in terms of helicity conserving and non-conserving 
amplitudes as follows: 

= (3/2) M 1+ + (1/2) E 1+ 
= E 1+ - Mi+ 

,1+ = (2Q 2 /pl)S 1+1 (44) 

where p* is the cm. pion momentum. Thus, helicity conservation implies B\ + = 0, or E±+ = M\ + . This is 
quite different from the low Q 2 situation. 
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QCD sum rule techniques were applied in [ 100 1 , to calculate the distribution functions for the A(1232) excita- 
tion. The CZ 1 32 proton wave function yields a small transition form factor, Q 4 Gt(Q 2 ) ~ 0.07 asymptotically. 
This can be traced to a cancellation in the leading order term of the matrix elements connecting the sym- 
metric A(1232) distribution amplitude, with the symmetric and antisymmetric proton distribution amplitude 
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Figure 11: The quantity Gt /Gdipoie verses Q 2 for the elastic form factor (a), and for transitions to the first (b), 
second (c) and third (d) resonances respectively, with Gdipoie = 2.79(1 + Q 2 /.71)~ 2 in (a) and 3(1 + Q 2 /.71)~ 2 
in (b-d). The first resonance (b) is the A(1232) (the Pa3(1232)). The second resonance ( c) at Q 2 above about 
3 GeV 2 is dominated by the Si2(1535). The third resonance at low Q 2 is dominated by the Fi 5 (1680). The fits 
for Gt were based on inclusive data referenced in 1 93 , 9^] . and selected data from |9{| . The elastic proton form 
factor Gmp is shown in a). Also shown at lower Q, denoted by (x), are form factors derived from amplitudes 
obtained from exclusive (e,e',p)n° and (e,e',p)rj data. The dashed curves are the result of the local duality 
sum rule calculations of and Q] for the elastic and A(1232) transitions respectively. The solid curve in a) 
is the Gmp result of the PQCD sum rule calcutation of |57j employing 4>cz- The solid lines at the lower right in 
b) and c) are the result of the PQCD calculation of |3^] using 4>cz for the A(1232) and 512(1535) respectively. 
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respectively. Schematically, | (0a \Th\4>s) + (0A|Tff|</>^)| is much smaller than either alone. 

If the leading amplitude of the P — > A transition is indeed small, the anomalous shape of the transition form 
factor might be explained as follows. At high Q 2 , the leading order helicity conserving amplitude dominates 
over the helicity non-conserving amplitude. That is, A 1 / 2 >> A 3 / 2 . A suppression of the A 1 / 2 amplitude at 
all Q 2 due to the cancellation of the symmetric and antisymmetric matrix elements, might then result in the 
dominance of the A 3 / 2 amplitude over a larger range of Q 2 than otherwise expected, and Q 4 Gt{Q 2 ) would 
decrease as a function of Q 2 . In fact the evidence that E\ + /Mi + is still small for Q 2 up to 3 GeV 2 , is consistent 
with the dominance of non-leading processes. 

Recently Q the local duality procedure was applied to the A(1232) form factor, and it was found, as in the 
pion case, that the form factor in the few GeV 2 region can be accounted for by purely soft processes (see Fig. 
|To| ) . However, it then falls significantly below the experimental values at higher Q 2 , which might be evidence 
that hard processes are playing an increasing role. 

It will be interesting in the future to determine whether Q 4 Gt(Q 2 ) does indeed level off above Q 2 = 10 
GeV 2 , and where the Ei + amplitude becomes comparable to the Mi+. This would support the valence PQCD 
description. 



3.3.3 The Second Resonance. Fig. 11 shows that at high Q the form factor for the peak at W ~ 1535 



MeV approaches the Q~ 4 dependence predicted by valence PQCD. Although the Ui3(1520) is dominant at Q 2 
= 0, there is a crossover and the 5n(1535) dominates the Di3(1520) at Q 2 ~ few GeV 2 |p7| . Another unique 
feature of the S , n(1535) is that it is the only excited state with a large n decay branching ratio (~ 50%), so 
that experimentally it is easily is isolated. 

Ref. |l00[ presents a calculation of the proton — > Sti transition form factor in the valence PQCD framework. 
The result is a behavior similar to the elastic form factor. Although the results are about a factor of two lower 
than the data, the authors remark that theoretical uncertainties in the distribution functions, and higher order 
contributions to a s are probably great enough to account for these discrepancies. 

3.3.4 The Third Resonance. Fig. |il"| shows that at high Q 2 the form factor for the peak near W = 1680 
MeV is consistent with the predicted Q behavior. The errors are large, however, and it is not clear how many 
resonances are contributing to this peak. The potential for obtaining separated resonance amplitudes at high 
Q 2 with exclusive reactions is very good. This is particularly true since it has been demonstrated [101, ^3 ; 102], 
that the non-resonant background diminishes with Q 2 at approximately the same rate as the resonances. 

3.3.5 Duality. A very interesting concept is that of duali ty b etween resonances and the non-resonance 
continuum in the W region where they overlap. One observes [ 101 1 that the rate of decrease with Q 2 of the 
resonance cross sections approximately follows the extrapolation of deep-inelastic scaling into the resonance 
region, suggesting that both processes are related by the same underlying physics. Later, this was put on 
firmer ground, and it was shown that leading logarithmic corrections extend the duality range in Q 2 [103, 104]. 
However, all of this is based on analyses of inclusive data, which cannot effectively separate non-resonance 
from resonance contributions. In order to access this very fundamental result one really needs to have a clean 
separation of resonance and non-resonance data over a large range of Q 2 , which can only be accomplished by 
the measurement of exclusive reactions. 



4 Conclusions 

As seen, much work remains in both experiment and theory. Valence PQCD and factorization appear to be an 
attractive starting point for treating high Q 2 form factors, although how high Q 2 must be for valence PQCD to 
dominate remains controvertial, and most probably depends on the specific reaction. Opinions on this matter 
vary strongly, from those who maintain that the required Q 2 is much higher than is likely to be experimentally 
accessible in the forseeable future, to those who believe that valence PQCD is already applicable at Q 2 as low as 
a few GeV 2 . We suggest that the quality and extent of existing data does not allow a definitive conclusion, but 
that soft non-perturbative processes probably play an important role for much of the existing data. On the other 
hand, given the complexity of QCD, there is a need for further theoretical work, based on fundamental principles 
of QCD, to deal with the soft, or Feynman mechanism. Indeed, it may be possible to express form factors in the 
transition region as a sum over valence PQCD and Feynman mechanism contributions, with the hard-scattering 
imbedded in the latter treated with PQCD methods. For the truly soft region, lattice calculations may play an 
increasing role in the future |2(| . 

On the experimental side, one must push the frontiers to as high Q 2 as technically feasible, to provide data 
which has the best chance of testing these ideas. One should also go beyond the experiments which merely 
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test constituent scaling, to those which test other central tenets of theory, such as helicity conservation. Such 



work has now begun at CEBAF, and may be further extended by a proposed European facility, ELFE [ 105 In 
summary, this area appears to offer some of the most interesting theoretical and experimental challanges for 
the next decade. 
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